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ABSTRACT 

In studies of the environmental dependence of structure formation, the large scale 
environment is often thought of as providing an effective background cosmology: e.g. 
the formation of structure in voids is expected to be just like that in a less dense 
universe with appropriately modified Hubble and cosmological constants. However, 
in the excursion set description of structure formation which is commonly used to 
model this effect, no explicit mention is made of the effective cosmology. Rather, this 
approach uses the spherical evolution model to compute an effective linear theory 
growth factor, which is then used to predict the growth and evolution of nonlinear 
structures. We show that these approaches are, in fact, equivalent: a consequence of 
Birkhoff 's theorem. We speculate that this equivalence will not survive in models where 
the gravitational force law is modified from an inverse square, potentially making the 
environmental dependence of clustering a good test of such models. 
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1 INTRODUCTION 

One of the standard predictions of nonlinear hierarchical 
structure f ormation models is the abundance of virialized 
structures JPress fc Schechter IOtI : ISheth fc Tormen 19991 : 
I Jenkins et al. 2001 ). Simulations show that this abundance 
depends on the large scale environment: the ratio of mas- 
sive to low mass objects is larger in dense regions (e.g., 
iFrenk et al. 19881 ). Recent measurements in galaxy surveys 
appear to bear this out: the virial radii of objects in un- 
derdense regions are smaller, consistent w ith their having 
smaller masses (e.g.. lAbbas fc Sheth 20071 ). 

This paper is motivated by the fact that there are 
currently in the literature three methods for estimating 
how the mass function of virialized halos depends on the 
environment which surrounds them. The first, and per- 
haps easis t to implement, is based on the excursi on set 
approach (|Mo fc White 19961 : ISheth fc Tormen 20021 ). The 
second argues that halos which form in, say, voids should 
be thought of as forming in a less dense background cos- 
mology, so the mass function is that in a univers e with 
f^void = ^^o(l + Avoid) (e.g., iGottlober et al. 20031 ). The 
third is similar, but notes that to correctly estimate the 
background cosmology, one must account not only for the 
lower density in a void, but for the fact the effective Hubble 
constant of t he void cosmology is large r than in the back- 
ground (e.g.. IColdberg fc Vogelev 2004 ). One way of think- 



E-mail: shethrk@physics.upenn.edu 



ing about the effective Hubble constant is that it ensures 
that the effective cosmology has the same age as the back- 
ground cosmology. (The cosmological constant is, of course, 
constant, but when expressed in units of the critical den- 
sity in the effective model, it is modified because the crit- 
ical density depends on the effective Hubble constant.) In 
Section (2] we use the spherical evolution model to show 
that the first a nd third methods are equivalent (although 
iGoldberg fc Vog elev 2004, state otherwise), and that both 
are incompatible with the second method (which incorrectly 
ignored the change to the Hubble constant). 



There has been recent interest in the fact that the for- 
mation hi stories of halos of fixed mass depend on their envi- 
ronment ISheth fc Tormen 20o3 : lGao et al. 20051 ). an effect 
which is not predicted by the simplest excursion set meth- 
ods (e.g.. White 199(|). So one might have wondered if this 
is where the difference between the excursion set approach 
and one based on the effective cosmology is manifest. In Sec- 
tion[3]we show that in this case also, the two approaches are 
equivalent. 



A final section summarizes our results, and speculates 
that the equivalence we have shown will not survive in mod- 
els models where the force law has been modified from an 
inverse square. 
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2 THE CALCULATION 

The main point of the foUowing calculation is to show ex- 
plicitly that, at least for cosmologies with no cosmologi- 
cal constant, the environmental dependence of halo abun- 
dance s can be described using the excursion set approach 
(e.g., IMo fc White 19961 : ISheth fc Tormen 2002l l. Namely, 
one need not worry about the details of the effective cosmol- 
ogy as sociated with the region sur rounding the perturbation 
(as do lGoldberg fc Vogelev 2004 ): it is enough to compute 
an effective growth factor using the spherical collapse model. 
Although we have phrased our discussion in terms of an 
flo = 1 background cosmology, it is obviously applicable to 
arbitrary values of no- Our analysis suggests that this re- 
mains true when the background cosmology has A 7^ 0. 

For what follows, it is useful to recall that the age- 
redshift relation in an Qo = 1 cosmology is given by 
H{z)t{z) = 2/3, where H is the Ifubble constant. In an 
open universe, this relation is 
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1, so a{t) = (1 + z)~^ , and 



with the co nvention that a (to 
Qo = ^{to) (|Peebles 1980[ ). 

The linear theory growth factor is D{t) 
and if f2o < 1 then 



a{t) if ^0 = 1, 



where x = a{t) (1 - D.o)/no (jPeebles lOSd ). 
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2.1 The spherical evolution model 

The spherical evolution model describes the evolution of the 
size i? of a spherical region in an expanding universe: 

d^R _ GM{< R) 

It provides a parametric relation between the 
sity contrast predicted by linear theory (5 (to) 
D{to)/D{tinit) S{tinit), the nonlinear overdensity A, 
the infall speeds Vr... (ICunn fc Gott 19721: [Schechter 1980l : 
[Peebles 1980l : iPadmanabhan 19931 : iBernardeau et al. 200^ . 
Here D{to) is the linear theory growth factor at time to, and 
we will often use the shorthand. So = S{to). 
If n = 1, then 
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where the first expression in each pair is for initially over- 
dense perturbations and the second is for underdense ones. 
Overdense perturbations eventually collapse, the final col- 
lapse being associated with the value 9 = 2n, at which time 
the hnear theory density is Sd = (3/5)(67r/4)2/3 = 1.68647. 
In this section, we use the subscript 1 to indicate that this 
value is associated with f2o = 1. 

If Slo < 1, then only perturbations above some density 
Smin will collapse, and 

f{9) t;pec _ 9(0) , 
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+ 1, 



arccos(2/r2o — 1) if closed 
arccosh(2/f2o — 1) if open ' 

is the Hubble constant, and 

9 sinh a; (sinhtj — uj) 
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Complete collapse is again associated with 9 = 2n, and we 
will write the critical linear density required for collapse as 



- Sd/H'^)]- 



(9) 



It happens that Scui depends only weakly on Q,o. When f2 ^ 
1, then Scu. (3/5)(67r/4)2/3 = 1.68647, and Scu. ^ 3/2 
when f2o 0. 

The parametric solution is rather cumbersome. It hap- 
pens that the relation between So and A is rather well ap- 
proximated by 



1 + A^ {l-So/5,^y 



(10) 



Similarly, it is also useful to have an approximation to the 
exact solution for the linear theory growth factor. When 
f2o < 1, then the linear theory growth factor is well approx- 
imated by 

(5/2) a(t) ^(t) 



(11) 



nity/"^ + i + Q.{t)/2' 

(jCarroll et al. 19921 ). where a(t) denotes the expansion fac- 
tor at time t, and Q.{t) is given by equation ([2]). This expres- 
sion is normalized so that D{to) = a (to) = 1 if ^0 = 1. 

2.2 Environment and spherical evolution 

Suppose we consider the evolution of a spherical underdense 
region in an fio = 1 universe. Let 1 -I- A^j < 1 denote the 
density in this region. If we wish to think of this region as 
being an underdense universe, then the effective value of 
Q in this region is smaller than unity for two reasons: first, 
because the density is lower by a factor of H- A^, , and second 
because the region is expanding faster than the background, 
so it has an effective Hubble constant H^i which is larger. 

To see what equation ((6)l implies for the evolution, let 
1 -I- Al denote the density of a small patch respect to the 
background density (the subscript unity denotes the fact 
that this is the overdensity with respect to a background 
which has critical density: f2o = 1). Now, suppose that this 
patch is surrounded by a region U within which the average 
density is l-f A„ with respect to the true background. Then 
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the smaller patch has overdensity (1 + Ai)/(1 + A^^) with 
respect to its local background. If we wish to describe the lo- 
cal environment as has having its own effective cosmological 
parameters, then the local value of the Hubble constant Hui 
differs from the global one Ho'- H^j/Ho ~ g{i^)- Thus, the 
expressions in equation © are really the statements that 



1 + Ai 



and 



"i^pecl Upecl 



(12) 



where ttpcci is the peculiar velocity of the shell U with 
respect to the background, had the mass within U been 
smoothly distributed (we know it is not because the cen- 
tral region has density 1 -I- Ai). Now, the local value of Q^j 
within U differs from the global value f2o = 1 both because 
A„ 7^ and because the different expansion rate means that 
the local value of the critical density is different: 



no(l + A„) _ f{uj) 



(13) 



where we have used the fact that Q.q = 1. Notice that this 
relation between 0,^^ and uj is the same as equation ((Tjl. In 
other words, we get the same description for the evolution 
of the small scale patch if we treat it as having overdensity 
1 + Ai with respect to the f^o = 1 background within which 
the Hubble constant is Hq, as if we describe it with respect 
to the local cosmological model Q.^^ and H^^, and we rescale 
our definitions of density and peculiar velocity accordingly. 
In addition, using the exact expression for the age of the 
universe given above, we can see that these definitions also 
guarantee that to is the same in the both the background 
and the local cosmological model. 

If we write the linear theory overdensity associated with 
1 -I- Ai^ as 



h{uj) 



(14) 



then 



5o = ^ \h{e) 



&ci — 3u 



The term in square brackets is simply the difference in linear 
theory values for the background cosmology. If we think of 
this as an effective linear theory overdensity in the effective 
cosmology, then the prefactor is the effective linear theory 
growth factor. It is straightforward to verify that, indeed. 



(15) 



where Di is the growth factor in the background cosmol- 
ogy, and Du, is the growth factor in the patch, at time to. 
This last point is important, as the expansion factor a{to) 
in the patch cosmology is not equal to the expansion factor 
in the background cosmology. In particular, we know that 
aui(to)/ai(to) = (1 + Auj)~^^^. For completeness, we note 
that 



1 + 



(2^) 



2/3 



(sinh L 



,)2/3 
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(recall that we are in an underdense region). 

In the following, take a-i{to) = 1, so Di = 1. The ap- 
proximate solution (|10|) of the spherical evolution model 



shows similar behaviour: 

1 + Ai n 



l + A 



Si/Sci 
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Scl — Su 



Suj/Scl 
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Scl — Si 

Scl — Sui 



Si-S^ 



(17) 



Scu^/D, 

where Si denotes the linear theory value associated with 
the nonlinear density Ai for fio = 1- The final approxi- 
mation follows from recalling that Scui depends only weakly 
on cosmology. Comparison with equation (|10p shows explic- 
itly that the relevant linear theory quantity is the difference 
between the S7o = 1 values for the perturbation and the 
environment, and this difference must be multiplied by the 
linear growth factor D^i in the effective cosmology. 

Now, to estimate the mass function of virialized ob- 
jects, we are interested in the case when 9 = 2n. The 
analysis above shows that 5cui/D^ = Sd — S^i; the ob- 
jects which form in a region of nonlinear density 1 -I- A^j 
with respect to the background, with corresponding lin- 
ear overdensity Sui, can either be thought of as forming i n 
an effective il^ cosmology (e.g.. IColdberg fc Vogelev 2004I '). 
or as forming in the true flo background cosmology but 
with an effective linear theory overdensity which is off- 
set by Sui to account for the surrounding overdensity (e.g., 
IMo fc White 19961 : [Sheth fc Tormen 20021 '). The second de- 
scription is easier to implement, and follows naturally from 
the excursion set description. In particular, the analysis 
above shows that approaches which do not correctly com- 
pute Q,uj (e.g., Gottlober et al. 2003 ignore the fact that 
Hui 7^ Ho) are incompatible with the excursion set ap- 
proach. In any case, the analysis above suggests that such 
approaches are ill-motivated. 



3 FORMATION HISTORIES 

The previous section showed that the excursion set approach 
results in the same expressions for the environmental depen- 
dence of the present day linear theory growth factor as one 
derives from thinking of the environment as defining an ef- 
fective cosmology. So the question arises as to whether or 
not the two approaches predict the same evolution. For ex- 
ample, one might have wondered if the formation histories 
of objects are the same in these two approaches. 

To see that they are, it will be convenient to modify our 
notation slightly. We showed that 

Sc{ilu.o) Sc{no) - Si^iAo) 
~n ^ n 

where the subscripts mean the present time. The quantity 
5l(Ao) is what we previously called 5^; it is the value of the 
initial overdensity extrapolated using linear theory (of the 
background cosmology) to the time at which the nonlinear 
density is Ao. Also, we previously had set the growth factor 
in the background universe at the present time to unity: 
Do = 1. We have written it explicitly here to show that, 
had we chosen to perform the calculation for some earlier 
time, then we would have found 



Sc{n^i) _ 5c(!^i)-5l(Ai) 



D^i 



Di 



(19) 
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where the subscript 1 denotes the earher time. I.e., Quii is 
the effective cosmology associated with the overdensity Ai, 
which itself is related to Ao by the spherical evolution model 
(the region that is Ao today was a different volume in the 
past, but its mass was the same.) And, analogously to the 
previous expression, (5l(Ai) is the initial overdensity extrap- 
olated using linear theory to the (earlier) time at which the 
nonlinear density was Ai. Since Ai is closer to than is Ai, 
5l(Ai) is also closer to than is 5l(Ao). 

If one were to apply the excursion set approach to study 
formation histories in the effective cosmology, one would be 
interested in the difference between equations (|19p and p8p : 



4(»i) _ 5,{no) 

Di Do 
"5l(Ai) 5l(Ao) 



Di 



Do 



(20) 



Now, the quantity in square brackets is 



5l(Ai) 5l(Ao) 



Di 



Do 



<5l(Ai 



Di/Do 



-5l(Ao) 



Do =0, (21) 



because (5l(Ai) and (5l(Ao) are the same quantity (the ini- 
tial overdensity) , evolved using linear theory to two different 
times. In particular, 5l(Ai) is closer to than is (5l(Ao) by 
5l(Ai)/5l(Ao) = Dr/Do. Thus, 



5c(nt^l) Sc{ftuio) 



D^ 



D^ 



Di/Do 



-5,{yio) 
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(22) 



Note that the expression on the right has no dependence on 
the effective cosmology. Moreover, it is exactly the same as 
the expression that one obtains when using the excursion 
set approach to study formation histories in the background 
cosmology. It is in this sense that the formation histories 
of objects are independent of the effective cosmology of the 
environment; the excursion set approach is a simple self- 
consistent way of exploiting this fact. 



4 DISCUSSION 

The excursion set description provides a simple, self- 
consistent way of estimating the effect of environment on 
structure formation and evolution. In particular, it is equiv- 
alent to using the fact that the large scale environment can 
be thought of as providing an effective background cosmol- 
ogy of the same age (Section [2]). Estimating the parameters 
of the effective cosmology is slightly more involved, but use- 
ful for running simulations which mimic the formation of 
structure in different environments. 

In essence, the equivalence between the excur- 
sion set and effective cosmology descriptions is a 
consequence of Birkhoff's theorem: the evolution of 
a perturbation does not depend on its surroundings. 
There has been recent interest in models with modi- 
fied gravitational for ce l aws (e.g.. IShirata et al. 20051 : 
IStabenau fc Jain 20061 : IShirata et al. 20071 ). Since 
Birkhoff' s theorem does not apply in such mod- 
els (iMartino et al. 20081: ISchafer fc Kovama 200"p : 



iDai et ah 20081 : IClifton 20061 : ICapozzieUo et al. 200'il ). 



it will be interesting to see if this equivalence survives. 
If not, the enviromental dependence of clustering may be 
added as another constraint on such models. 
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